Dynamics of quasi-periodic, bifurcation, sensitivity and three-wave solutions for (n + 1)-dimensional generalized Kadomtsev-Petviashvili equation

This study endeavors to examine the dynamics of the generalized Kadomtsev-Petviashvili (gKP) equation in (n + 1) dimensions. Based on the comprehensive three-wave methodology and the Hirota’s bilinear technique, the gKP equation is meticulously examined. By means of symbolic computation, a number of three-wave solutions are derived. Applying the Lie symmetry approach to the governing equation enables the determination of symmetry reduction, which aids in the reduction of the dimensionality of the said equation. Using symmetry reduction, we obtain the second order differential equation. By means of applying symmetry reduction, the second order differential equation is derived. The second order differential equation undergoes Galilean transformation to obtain a system of first order differential equations. The present study presents an analysis of bifurcation and sensitivity for a given dynamical system. Additionally, when an external force impacts the underlying dynamic system, its behavior resembles quasi-periodic phenomena. The presence of quasi-periodic patterns are identified using chaos detecting tools. These findings represent a novel contribution to the studied equation and significantly advance our understanding of dynamics in nonlinear wave models.


Introduction
It is widely acknowledged within the scientific community that nonlinear evolution equations (NLEEs) play a crucial role in the investigation of nonlinear physical phenomena [1][2][3][4][5].Contemporary research on multi-component integrable models has emerged as a significant advancement in the exploration of soliton equations [6][7][8][9][10].In recent years, there has been a proliferation of powerful techniques for investigating the exact solutions of the NLEEs through the use of symbolic computation [11][12][13][14][15].These methods include the modified simple equation method [16], the improved Bernoulli sub-equation method [17], Generalized logistic equation method [18], the Hirota direct approach [19,20], multivariate bilinear neural network method [21], variable-coefficient symbolic computation approach [22], and the Lie symmetry approach [23], among others.The Hirota bilinear method, which was developed by Hirota, stands out as a powerful and effective approach for constructing exact solutions of NLEEs [24].Once the bilinear form of a nonlinear equation is obtained by a dependent variable transformation, it becomes relatively effortless to obtain its multi-soliton solutions [24][25][26].Furthermore, by utilizing the Lie symmetry analysis and dynamical system method, one can also derive the symmetries and exact explicit solutions of NLEEs [27,28].
Complex system theory can be expounded as a mathematical branch that is devoted to scrutinizing dynamical systems that are characterized by a substantial number of variables [29].The fundamental aim of this field is to extend the principles of dynamical system theory, which primarily concerns systems that are composed of a few variables.They have been widely used in various fields such as robotics, engineering structures, and science [30][31][32].One illustrative example of a dynamical system is a framework utilized to encode trajectories from human demonstration, allowing for adaptability to changing environments and robustness to perturbations.Furthermore, the dynamical system of the hybrid electric vehicle incorporates front and back shaft drives, a battery, and a monitoring module which form an intricate system that requires the application of complex system theory to comprehend its dynamical behavior.
In the present study, we aim to provide a comprehensive and in-depth analysis pertaining to the behavior of both the simple and perturbed dynamical systems of the nonlinear physical model that represents a physical phenomenon.Our discussion encompasses a range of topics such as the bifurcation analysis, sensitivity analysis, and the irregular and unpredictable behavior that may be observed within the desired dynamical system.Additionally, we shall employ a variety of efficient techniques and tools to facilitate our investigation.
Bifurcation refers to a topological and qualitative alteration that occurs in the phase space of a system when certain parameters undergo a slight modification beyond their critical thresholds [29].Bifurcations assume significant functions in various tangible systems as a mode of switching.Such instances comprise the stimulation of neurons, the generation of patterns in morphogenesis (which will be expounded subsequently), the catastrophic shift of ecosystem states, and the storage of binary information in computer memory, among others.Seshasayanan et al. [33] investigate the primary bifurcations of a plane parallel flow with Kolmogorov forcing, including a new type called stationary drift bifurcation.Tang [34] delves into the exploration of bifurcation and dispersive soliton solutions within optical fibers utilizing the Schro ¨dinger-Hirota equation.
Chaos refers to the persistent behavior of a nonlinear dynamic system that never follows fixed or repeating paths over time [29].It may appear as a random variation, but it does occur within entirely deterministic and simple dynamic systems.Despite being subject to deterministic rules, chaotic systems manifest a remarkable susceptibility to initial conditions, rendering long-term prognoses and exact outcomes exceptionally arduous.In continuous time systems, chaos may arise in systems featuring no less than three independent dynamical variables, while simultaneously being nonlinear.These conditions must be met for the detection of chaos in continuous time systems.Whenever a system exhibits sensitivity to initial conditions, it is deemed chaotic.In essence, the definitive description of chaos states that a chaotic system is a deterministic one that displays behavior that appears random and unpredictable.The vast expanse of literature presents a plethora of methodologies for the identification of chaos [35].These methodologies include time series analysis, three-dimensional phase portraits, Poincare ḿaps, power spectrum analysis, bifurcation diagrams, and Lyapunov exponents.These tools have demonstrated significant effectiveness in distinguishing the chaotic, quasi-periodic, and periodic patterns exhibited by physical models.
The Kadomtsev-Petviashvili (KP) equation, which was introduced by Kadomtsev and Petviashvili in 1970 as an analogue of the classic Korteweg-de Vries equation in two spatial dimensions, exhibits a high degree of efficacy in modeling nonlinear phenomena in various fields including fluid physics, plasma physics, Bose-Einstein condensates, optics, and beyond.The KP equation has garnered significant scholarly interest due to its physical and mathematical significance.Through the utilization of singularity manifold analysis, both Dorizzi et al. [36] and Wazwaz et al. [37] have demonstrated the integrability of the KP equation in instances of both constant and variable coefficients.Tian and Gao have successfully derived a variety of families of travelling wave solutions of the KP equation through their research efforts, as documented in their publication [38].The lump solution, the lump and one stripe soliton, and resonance stripe soliton solutions for the KP equation were then subsequently constructed by utilizing the pioneering Hirota's bilinear method [39].Several research advancements have been accomplished in regards to the (3+1)-dimensional KP equation and its variants [5,19,[40][41][42][43][44][45][46][47][48].The KP system is a nonlinear partial differential equation that generalizes the Korteweg-de Vries (KdV) equation to two spatial dimensions.Its significance lies in its ability to model wave propagation across different physical domains, including shallow water waves and plasma physics.An in-depth comprehension of the outcomes linked to the KP system necessitates an analysis of solutions, stability, and the mathematical characteristics inherent in the system.This system represents a fertile research area within the realms of applied mathematics and theoretical physics, providing profound insights into the dynamics of nonlinear waves.
The proposed extended form of the KP equation in higher dimensions is presented as follows which is the generalized (n + 1)-dimensional KP equation.Ma et al. [49] have expanded on the KP equation in (n + 1)-dimensions by introducing a series of lump solutions produced from quadratic functions.This was followed by Chen et al. [50] who developed the multi-lump or lump-type solutions.Nevertheless, the extended KP equation is not integrable for n � 3.If we take x 1 = x, x 2 = y, α = 6, β = 1, γ = 3ρ 2 , σ 1 = σ 2 = 0, then Eq (1) reads as when the values of ρ are equal to i and 1, Eq (2) corresponds precisely to the KPI and KPII equations, respectively.The KP equation can be used to model water waves of long wavelength with weakly nonlinear restoring forces and frequency dispersion.The alteration in the sign of ρ 2 is linked to the respective values of the gravitational and surface tension forces.When the values of x 1 , x 2 , and x 3 are assigned to x, y, and z respectively, Eq (2) can be simplified to the generalized KP equation in a (3+1)-dimensional space.
It should be noted that a (3+1)-dimensional KP equation that was examined in [40] did not pass the test for integrability, as a result of the presence of a second-order dispersion term U zz .In 2016, Wazwaz and El-Tantawy [51] discover the multiple solitons of model using the simplified Hirota's direct method.They prove that the new model fails the Painleve integrability test although it gives multiple-soliton solutions.In 2017, Sonmezoglu et al. [52] investigated the exact solitary wave solutions of model.Also, Xu and Wazwaz [53] undertook a research endeavor on Eq (3) with the objective of investigating its integrability and localized solutions.Through the application of singularity manifold analysis and binary Bell polynomial method, it has been determined that the generalized KP equation of (n + 1)-dimensions contains N-soliton solutions, alongside the Painleve ´property, Lax pair, Ba ¨cklund transformation, and infinite conservation laws, thus establishing its complete integrability.The aim of this research endeavor is to fabricate the three wave interactions and employ the Lie symmetry approach in order to identify the symmetry reductions, bifurcation analysis, sensitivity analysis and chaotic phenomena that arise in the intended dynamical system for Eq (3).
The subsequent sections of the paper are organized in the following order.In Section (2), the Hirota bilinear form that corresponds to the governing equation is presented.Section (3) serves as a tool for the construction of three-wave interactions, which are derived from the bilinear form presented in Section (2).In Section ( 4) and ( 5), we determine the symmetry reductions of the studied equation to obtain second order ordinary differential equation (ODE) using the Lie symmetry approach.Sec. ( 6) and ( 7) are used to study the bifurcation analysis and quasi-periodic behavior of the dynamical system.In Sec. ( 8) is used to perform sensitivity analysis using three different initial conditions.The main findings of the study are discussed in detail in Sec.(9).Sec.(10) provides the concluding remarks of study.

Hirota bilinear form of the governing equation
Consider the transformation defined as Uðx; y; z; tÞ ¼ 12b a ðln gÞ xx : ð4Þ Switching Eq (4) into Eq (3), we arrive at the following bilinear form as where the Hirota bilinear operator is defined as Using the above formulas of Hirota operator, we acquire 3 Three-wave solutions The following test function is used to ensure that the solutions to Eq (3) has three-wave solutions where ρ i = j i x + k i y + l i z + m i t, i = 1, 2, 3 and j i , k i , l i and m i are constants to be determined.By inserting Eq (8) for Eq (7) and equating like coefficients of different powers of e r 1 ; e À r 1 ; cosðr 2 Þ; sinðr 2 Þ; cosðr 3 Þsinðr 3 Þ and constant term to zero, we get: We can obtain three-wave solution of the (n + 1)-dimensional generalized KP equation by using Eq (8) with above known parameter values as follows: 4 Lie symmetry analysis of Eq (3) In this part, the aim is to identify the symmetry transformations of the studied equation.This is achieved by employing the well-known technique called Lie symmetry analysis, which was introduced by Marius Sophus Lie.Lie's revolutionary ideas transformed the exploration of continuous symmetries and offered a potent approach to comprehending and resolving differential equations.Nowadays, Lie symmetry analysis persists as a dynamic field of study, discovering novel applications across diverse disciplines.Let us examine the Lie group of infinitesimal transformations in (x, y, z, t, q), which is characterized by a single parameter.
with group parameter ε � 1 and ρ 1 , ρ 2 , ρ 3 , τ and z are coefficient functions.The vector field which is associated with the previously mentioned Lie group can be described in the following way: V ¼ r 1 ðx; y; z; t; UÞ @ @x þ r 2 ðx; y; z; t; UÞ @ @y þ r 3 ðx; y; z; t; UÞ @ @z þtðx; y; z; t; UÞ @ @t þ zðx; y; z; t; UÞ @ @q ; ð10Þ where the coefficient functions r 1 ðx; y; z; t; UÞ, r 2 ðx; y; z; t; UÞ, r 3 ðx; y; z; t; UÞ, tðx; y; z; t; UÞ and zðx; y; z; t; UÞ are to be evaluated later.The invariance condition for Eq (4) with V becomes: where Pr ð4Þ V is a fourth prolongation of V and is defined as: By employing the prolongation defined in Eq (12) and conducting a coefficient comparison for the different derivatives of the dependent variable, we get a linear system of PDEs that is over-determined.Subsequently, after carrying out some computations, we get the set of determining PDEs.After solving the system of PDEs, we are able to successfully obtain the following set of symmetry transformations.
where F 1 (t), F 2 (t) and F 3 (t) are free arbitrary functions.Symmetry transformations given in Eq (13) can be used to reduce the dimension of the considered equation.Therefore, in the next portion, we determine the reduced dimensions of the governing equation using computed symmetry transformations.

Symmetry reductions
As mentioned in the previous section, we have to determine the reduced dimensions of the governing equation using the symmetry transformations given in Eq (13).The reason behind to introduce this section is that we can obtain an ODE which is further used to investigate the qualitative analysis and chaotic structures of the considered equation.Reduction using subalgebra: We assume the characteristic equation for V 1 as The desired solution for the characteristic Eq ( 14) is Using the Eq (15) in Eq (3), we get the following result: where F ðx 2 ; x 3 Þ and Gðx 2 ; x 3 Þ are arbitrary functions of ξ 2 and ξ 3 .Reduction using subalgebra: We assume the characteristic equation for V 2 as The desired solution for the characteristic Eq (17) is Using the Eq (18) in Eq (3), we get the following PDE: Eq (19) has the following three translation symmetries: ( Symmetry η 1 + η 2 + η 3 yields the following invariants: Substituting Eq (20) into Eq (19), we have The infinitesimal generators of Eq (21) are with c 1 , c 2 , c 3 and c 4 arbitrary constants.Therefore, the Lagrange system of Eq ( 22) is The group invariant solution of Eq ( 23) for c 1 6 ¼ 0 and c 2 = c 3 = c 4 = 0 is Inserting Eq (24) into Eq (21), we get Reduction using subalgebra: The desired solution for the characteristic Eq (26) is Using the Eq (27) in Eq (3), we get the following PDE: Eq (28) has the following three translation symmetries: ( Symmetry s 1 η 1 + s 2 η 2 + η 3 yields the following invariants: Substituting Eq (29) into Eq (28), we have The infinitesimal generators of Eq (30) are with c 1 , c 2 , c 3 and c 4 arbitrary constants.Therefore, the Lagrange system of Eq (31) is The group invariant solution of Eq (32) is Inserting Eq (33) into Eq (30), we get For convenience, we take γ = 0 in the above equation and it results into Eq ( 35) is seen to be of significant importance and has the potential to be used in further studies.Now, using the Galilean transformation to Eq (35), we obtain the subsequent dynamical system where

Hamiltonian saddle-node bifurcation
Hamiltonian saddle-node bifurcations offer a key mechanism for exploring how stationary points evolve within a Hamiltonian system in response to parameter changes.Hamiltonian saddle-node bifurcations are of paramount importance in comprehending and forecasting the dynamics of various physical and engineering systems.These bifurcations are instrumental in elucidating the evolution of systems with varying parameters, providing essential insights into the stability, transitions and qualitative dynamics of Hamiltonian systems.For this, the Hamiltonian function corresponding to system (36) is given as: where h denotes the Hamiltonian constant.To find the equilibrium points of the dynamical system (36), we need to solve the following system of equations: The equilibrium points of the system (36) are: Furthermore, the Jacobian matrix of the linearized system ( 36) is computed and its determinant is given as As we know Based on the above points, the subsequent observations can be deduced concerning the real axis line ϕ: a. (ϕ, 0) is a saddle point when The various outcomes that arise from varying the involved parameters are as follows: Case 1: For d 1 > 0 and d 2 > 0 By choosing the suitable parameter values λ 1 = 1, λ 2 = β = 0.5 and With reference to the present set of parameters, it can be observed that the two points of equilibrium are situated at coordinates (0, 0) and (1, 0), respectively.The visual depiction of their respective phase portraits can be found in Fig 2(a).These portraits show that (0, 0) is a saddle point, whereas (1, 0) is a centre point.
Case 2: For d 1 < 0 and d 2 > 0 By choosing the suitable parameter values λ 1 = −0.5, λ 2 = β = 0.5 and With reference to the present set of parameters, it can be observed that the two points of equilibrium are situated at coordinates (0, 0) and (−2, 0), respectively.The visual depiction of their respective phase portraits can be found in Fig 2(b).These portraits show that (0, 0) is a center point, whereas (−2, 0) is a saddle point.
Case 3: For d 1 > 0 and d 2 < 0 By choosing the suitable parameter values With reference to the present set of parameters, it can be observed that the two points of equilibrium are situated at coordinates (0, 0) and (−1, 0), respectively.The visual depiction of their respective phase portraits can be found in Fig 2(c).These portraits show that (0, 0) is a saddle point, whereas (−1, 0) is a centre point.
Case 4: For d 1 < 0 and d 2 < 0 By choosing the suitable parameter values λ 1 = −0.5, With reference to the present set of parameters, it can be observed that the two points of equilibrium are situated at coordinates (0, 0) and (2, 0), respectively.The visual depiction of their respective phase portraits can be found in Fig 2(d).These portraits show that (0, 0) is a center point, whereas (2, 0) is a saddle point.

Quasi-periodic behavior of perturbed autonomous system
To investigate the quasi-periodic and chaotic dynamics of the governing equation, we introduce a perturbation term cos(ωs) into the dynamical system (36).Consequently, the perturbed autonomous dynamical system of (36) with external periodic force can be written as follows: where θ = ωs.In system (40), the parameters denoted by g 0 and ω correspond to the frequency and magnitude of the external force applied, respectively.The phenomenon of chaotic solutions in dynamical systems is a captivating subject that has received thorough investigation within the realm of nonlinear dynamics.The term "chaos" denotes a behavior distinguished by its sensitive reliance on initial conditions, its aperiodic and unpredictable trajectories, and the presence of a strange attractor within phase space.The identification and verification of chaotic structures in the perturbed dynamical system has been executed by means of chaos detecting tools including three dimensional phase portrait, Poincare map, time series analysis, Lyapunov exponents, as discussed in introduction section.

Sensitivity analysis
The importance of conducting sensitivity analysis on a dynamic system is rooted in its ability to provide valuable insights into how slight modifications in initial conditions or parameters can impact the behavior of the system.Such analysis is crucial in understanding the stability and predictability of intricate systems, and can assist in making informed decisions and formulating effective control strategies.In this section, we undertake a comprehensive study of the sensitivity of a dynamical system (36) at three distinct initial conditions, as illustrated in Fig 7.This analysis is performed under the suitable choice of parameter values and these are d 1 = 2 and d 2 = 3.5.It is observed that the dynamical system (36) is not highly sensitive; small variations in its initial conditions do not result in significant deviations from the path followed by the solution.

Results and discussion
The current study involves development of new systems that explore useful findings in solitary waves theory.In Sec.
(3), we discover the three-wave solutions using anstaz function and present them in 3D and density plots as shown in Fig 1 .Sec. ( 6), signifies the dynamical observations of the planar dynamical system using bifurcation analysis.We show that how this dynamical system can be affected by its parameter values within the system, see Fig 2.
Further, we introduce the perturbation term to show the quasi-periodic behavior in the nonlinear dynamical system using different chaos detecting tools available in the literature as shown in Figs 3-7.In Fig 3, the periodic and quasi-periodic behavior of the solution is observed using phase portrait tool.In Fig 4, we use Poincare section tool to identify the periodic and quasi-periodic behavior using g 0 = 1 and g 0 = 1.5.In Fig 5, time series analysis is performed for perturbed dynamical system at time span [0, 100] to observe the quasi-periodic behavior of solution.In Fig 6, we calculate the Lyapunov exponent values for a perturbed dynamical system at g 0 = 0.5 and g 0 = 2.5, given in Table 1.It also includes the sensitivity analysis of the planar dynamical system (36) at different initial values and observe the behavior of solutions as can be seen in Fig 7 .The derived solutions could be very significant in elaborating physical aspects of real-world phenomena.The consequences of bifurcation and chaotic dynamics in dynamical systems encompass significant alterations in the qualitative behavior of the system, whereby minor perturbations in initial conditions can give rise to markedly diverse outcomes, rendering long-term forecasting challenging.Bifurcation may lead to the abrupt emergence of novel, stable states, whereas chaos introduces unreliability and susceptibility to initial conditions.These phenomena can have a substantial influence on domains such as meteorology, engineering, and financial markets, where grasping and handling the inherent unpredictability is essential for modeling and decisionmaking.

Conclusion
The work shows that the standard KP equation, and a variety of its extensions, give novel of integrable and non-integrable systems that lead to distinct physical structures.The work shows that this widely used model, that describes scientific features in plasma physics and fluid dynamics, explores many physical properties when extended to other versions.Firstly, we have applied the extended three-wave approach and Hirota's bilinear method to obtain novel threewave solutions.These solutions have been portrayed through figures to demonstrate an abundance of physical structures, as shown in Fig 1 .Next, the Lie symmetry approach has been applied in order to elucidate the symmetry transformations that facilitate the reduction of equation dimensionality.Through utilization of the subalgebra of these transformations, a second-order differential equation has been successfully obtained.After that, we apply the Gallian transformation to the derived ordinary differential equation given in Eq (35).Sec. ( 6) offers an analysis of Hamiltonian saddle-node bifurcation and sensitivity of the dynamical system (36).
In Fig 2, we have presented the phase portraits of Hamiltonian saddle-node bifurcation for different case of d 1 and d 2 .To observe the quasi-periodic motion of the governing equation, we consider the perturbed dynamical system (40) involving perturbation term.When an external force is applied to the system (36), its behavior displays characteristics that bear similarity to those observed in chaotic phenomena.The identification of chaotic behavior in perturbed dynamical system is analyzed using phase portraits, Poincare sections, time series analysis and Lyapunov exponents as shown in Figs 3-6.Sensitivity analysis of the governing equation has done and it is observed that (see Fig 7) small variations in initial conditions of the dynamical system (36) do not result in significant deviations from the path followed by the solution.The derived solutions could be very significant in elaborating physical aspects of real-world phenomena.The final remarks presented in this article provide 325 additional insight into the evolutionary progression and physical mechanisms governing 326 solitons, thereby augmenting our comprehension of the nonlinear dynamics 327 characterizing solitons in optical communication systems.